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Abstract. We present two-dimensional (2D) radiative transfer modeling of IRC + 10 216 at selected mo ments of its 



evolution in 1995-2001, which correspond to three epochs of our series of 8 near- infrared speckle images (Osterbart 



.it al. 2000 



Weigelt et al. 2002). The high-resolution images obtained over the last 5.4 years revealed the dynamic 



evolution of the subarcsecond dusty environment of IRC +10 216 and our recent time-indep endent 2D radiative 



trans fer modeling reconstructed its physical properties at the single epoch of January 1997 (Men'shchikov et al 



2001). Having documented the complex changes in the innermost bipolar shell of the carbon star, we incorporate 
the evolutionary constraints into our new modeling to understand the physical reasons for the observed changes. 
The new calculations show that our previous static model is consistent with the brightness variations seen in the 
near- infrared images, implying that during the last 50 years, we have been witnessing an episode of a steadily 
increasing mass loss from the central star, from M « 10 -5 Mq yr _1 to the rate of M ~ 3 x 10 -4 Mq yr _1 
in 2001. The rapid increase of the mass loss of IRC +10 216 and continuing time-dependent dust formation and 
destruction caused the observed displacement of the initially faint components C and D and of the bright cavity 
A from the star which has almost disappeared in our images in 2001. Increasing dust optical depths are causing 
strong backwarming that leads to higher temperatures in the dust formation zone, displacing the latter outward 
with a velocity vt ~ 27 kins" 1 due to the evaporation of the recently formed dust grains. This self-regulating 
shift of the dust density peak in the bipolar shell mimics a rapid radial expansion, whereas the actual outflow 
has probably a lower speed v < Voo ~ 15 kms -1 . The model predicts that the star will remain obscured until M 
starts to drop back to lower values in the dust formation zone; in a few years from that moment, we could be 
witnessing the star reappearing. 
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1. Introduction 



Near-infrared images with resolutions better than 100 



mas presented by (Weigelt et al. 1997; |Osterbart et al. 


The prjlsating carbon star IRC +10 216 (also known as CW 1997 


; Weigelt et al. 1998a 




Haniff & Buscher 1998; 


Leo, AFGL 1381), together with its huge circumstellar en- psterbart et al. 2000|; Tuthill et al. 2000; 


Weigelt et al. 



velope lost during its long evolution on the asymptotic gi- 
ant branch (AGB), is the best studied object of its class. 
Being in a very advanced phase of its life, probably in tran- 
sition to protoplanetary nebulae ( Osterbart et al. 2000| , 
hereafter Paper I), IRC +10 216 presently exhibits a very 
high mass-loss rate M ~ 10 -4 MQyr -1 (Men'shchikov 
et al. 2001, hereafter Paper II). After three decades of 



intensive o bservatio nal and theoretical work (see, e.g., ref- 
erences in Paper II] ), recent near-infrared speckle imaging 
has revealed an extremely complex evolution of its circum- 
stellar material in the vicinity of the dust condensation 
zone on a time scale of one year. 



Send offprint requests to: A. B. Men'shchikov 



2002) have demonstrated that the inner, subarcsecond 
dust shell of IRC +10 216 is non-spherical and clumpy, 
with four components A, B, C, and D clearly visible. The 
detailed two-dimensional radiative transfer modeling pre- 
sented in Paper II has shown that the star is actually 
located at the position of the second brightest component 
B. The brightest southern peak A was identified with the 
radiation emitted and scattered in the optically thinner 
cavity of the dense circumstellar shell. The model recon- 
structed physical properties of the star and dusty enve- 
lope of IRC +10216 at a single moment corresponding 
to the epoch of our high-re solution H- and if-band im- 
ages on January 23, 1997 ( |Paper l| ). As much as it was 
possible, the model took into account most other obser- 
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vations of dust radiation, although stellar pulsations and 
non-periodic changes of the shell made many measure- 
ments from various epochs fundamentally incomparable 
in the frame of the static model. 

In the present study, we attempt to attack the problem 
using a simplified approach based on the self-consistent 
model fo r a single epoch that we have constructed in 
Paper II . The idea is to extend this modeling to the first 
and to the last epochs of the 6-year sequence of high- 
resolution if-b and speckle images w e have obtained since 
October 1995 flWeigelt et al. 2QQ2| , hereafter Paper III). 
In the beginning of our monitoring of IRC +10 216, com- 
ponent B (the star) was relatively bright. The direct 
light from the star has been gradually fading since then, 
whereas the bright lobe A (the southern cavity) has be- 
come dominant. The angular distance between compo- 
nents A and B increased from ~ 190 to ~ 350 mas between 
1995 and 2001, implying a linear speed va ~ 18kms _1 
in the plane of sky and a deprojected radial velocity 
v r ,A ~ 19kms _1 (Appendix A of Paper II), for the as- 
sumed distance of D = 130 pc and the viewing angle of 



V = 40° ( [Paper lip . On the basis of our time-independent 
model, this well-documented rapid evolution has been 
qualitatively interpreted in terms of dust formation in a 
bipolar stellar wind, the increasing mass-loss rate, and the 
sublimation of the recently formed grains in a progres- 
sively hott er region just outside the dust formation radius 
([Paper III|) . 

The goal of the present study is to determine, whether 
the observed changes in the high-resolution speckle images 
are consistent with our previous model of IRC +10 216, 
and to derive more accurate physical parameters of the 
wind at several moments in time. In Sect. ||, we describe 
our approach and model parameters. In Sect. ||, we discuss 
the results of our new modeling, comparing them with 
the high-resolution images, intensity profiles, and spectral 
energy distribution (SED) of IRC +10 216. In Sect. |, we 
summarize the model parameters and our conclusions. 



2. Radiative transfer model 

2.1. General formulation 

Our new 2D radiative transfer modeling is based on the 
structure and physical properties of IRC +10 216 derived 



in Paper II . We refer to the latter for a detailed description 



of the model geometry and discussion of all assumptions. 
The parameters of the envelope were kept unchanged, ex- 
cept for the density distribution within 100 AU from the 



star and the cavity's opening angle (see Sects. 2.2 and 3.1). 
The density profile of the circumstellar shell on subarcsec- 
ond scales was allowed to vary due to expansion, increased 
mass-loss rate, and condensation of new dust in the wind 



(Sects. |3jJ, I3J, and|3J). 

The new models were computed for selected epochs 
(1995, 1997, and 2001) spanning the entire sequence of our 
speckle images in the K band (Paper III). Slightly adjust- 
ing the density distribution of the time-independent model 



of Paper II , we found the models which would produce 
the images and intensity profiles consistent with the ob- 
served ones. We limited ourselves to the constraints from 
the high-resolution near-infrared images and did not at- 
tempt to perform a new global search of the entire param- 
eter space for each epoch in order to fit all observations. 
This approach is justified by the fact that the latter have 
already been explained self-consistently by our previous 
model. Moreover, it is highly problematic to use a mix- 
ture of old observational constraints in a meaningful way 
for the rapidly evolving inner shell of IRC +10 216. 

2.2. Model parameters 

To facilitate the presentation and understanding of the 
new modeling, we briefly describe the parameters of our 
previous, time-independent model of IRC +10 216. At the 
end of this section, we summarize the time-dependent pa- 
rameters of the new models. 



The study reported in |Paper II] reconstructed the fol- 
lowing physical picture. For the adopted distance D = 130 
pc, the luminosity L* of the pulsating star varies from 
13000 to 5200 Lq between the maximum and minimum 
brightness, with a period of 649 days. The effective tem- 
perature T* of the star is changing from 2800 to 2500 
K between the two phases, whereas the stellar radius R* 
varies in a relatively narrow range, from 500 to 390 Rq. 
During the last 1000 years, the star has experienced at 
least two episodes of very high mass loss with rates reach- 
ing M ~ 10 -4 Mq yr _1 . The broken-power-law density 
profile of our model (see Fig. ||; also Table 6 in Paper II ) 
results most likely from a complex interplay between stel- 
lar pulsations, changing mass- loss rate, and radiative ac- 
celeration of dust and gas in the innermost parts of the 
envelope. The latter has bipolar cavities with a full open- 
ing angle uo = 36°, tilted toward us by V = 40° from the 
sky plane. The huge envelope of a radius R2 = 3 pc has 
a total mass M = 3 Mq, which corresponds to a dust- 
to-gas mass ratio pa/p = 0.0039. The mass is consistent 
with an initial mass of the central star M® w 4 Mq and a 
present-day core mass M* ^0.7 Mq. Different dust com- 
ponents exist in the envelope of IRC +10 216 everywhere 
from the stellar photosphere to the outer boundary, pro- 
ducing the total visual optical depth ry « 40 toward the 
star, whereas in the polar regions it is r p y ^10. 

Both the radial distribution and composition of dust in 
the model of IRC +10 216 were strongly inhomogeneous. 
The innermost dust formation zone of the envelope was 
located very close to the stellar photosphere, where the 
newly formed SiC grains had temperatures T$ic ~ 2000 
K. Most of the carbonaceous dust nucleated further away 
from the star (at tq — 30 AU at that epoch), whereas 
the most volatile components condensed at even larger 
distances (at r[ Mg5Fe ]s = 45 AU). The abundances and 
size distribution of dust grains in our model depend on 
the distance from the star. The modeling identified the 
carrier of the 11.3 pm feature with [SiC, C] particles, i.e. 
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Fig. 1. Comparison of our new K-band model images (lower panels) with our K-band speckle images of IRC +10 216 
(upper panels) obtained in October 1995, January 1997, and March 2001 with resolutions of 92, 87, and 73 mas, 
respectively ( Paper III ). The two dashed circular segments show the radius of the density peak in the main dust 
formation zone of our model (Fig. ^) and the small dashed circles indicate the position and actual size of the star. 
The straight dashed lines intersecting at the star's center outline the conical bipolar cavities of the model. Each image 
shows a 1" x 1" area. North is up and east is to the left. 



unorganized aggregates made of incompletely amorphous 
carbon grains with significant graphitic content and silicon 
carbide grains. The 27 /im emission band, also known as 
the 30 fim feature, was well reproduced by the core- mantle 
[SiC,C] + [Mg,Fe]S grains. 

The new models of IRC +10 216 presented in this pa- 
per correspond to the epochs of October 8, 1995, January 
23, 1997, and March 9, 2001 (below, we refer to the epochs 
only by the years). We fixed the model parameters at 
the values derived in our previous model except for only 
the s tar's prop erties and the envelope's density distribu- 
tion ( [Paper l| ). The latter resulted from the most recent 
episode of increased mass loss (within a radius of ~ 100 
AU, Sect. |3^3|) . Thus, our model assumes that at the first 
epoch of our imaging series, in October 1995, the star had 
L*i = 12200 L , T*i = 2760 K, and = 483 R Q (17.3 
mas); at the third epoch, in January 1997, L*3 = 6140 
L , T* 3 = 2530 K, and R* 3 = 408 Rq (14.6 mas); and 
at the eighth epoch, in March 2001, L* 8 = 12900 L , 
T* 8 = 2790 K, and i?* 8 = 486 Rq (17 A mas). The cav- 
ities' opening angle changed from uo\ = 46° to uos = 36° 
to uo$ = 26° between the epochs (Sect. 3.1.3[ ). The latter 
correspond to the luminosity phases <\>\ = 0.88, ^3 = 0.61, 
and 08 = 0.93, respectively. The intervals between the 
epochs amount to Ati3 = 473 days and A^g = 1506 days. 

3. Results and discussion 



3.1. Model images and intensity profiles 

Figure [l| compares our new if -band model images with the 
high-resolution speckle images of IRC +10 216 obtained in 
1995, 1997, and 2001. Figure | displays the H-K, J-H, 
and J—K model color maps and the observed color distri- 
butions derived from the J, H, and K speckle images of 
2001. Figure || shows a more quantitative comparison of 
the model and observed intensity distributions in terms of 
the normalized intensity profiles. The latter were derived 
from the images in two directions, parallel (PA « 20°) 
and perpendicular (PA « 110°) to the projected symme- 
try axis of the bipolar structure. We define the axis as a 
straight line through the peaks A and B in the K image of 
1995, where the peaks are best visible. Before extracting 
the profiles, we convolved the model intensity distributions 
with the point-spread functions (PSF) of 93, 87, and 73 
mas full width at half-maximum (FWHM) corresponding 
to the resolutions of the observed images. 

To facilitate direct comparisons with observations, an 
outline of the model geometry is displayed in Figs. |l| and 
^| by dashed lines. The central star is shown by a small 
dashed circle whose diameter corresponds to the actual an- 
gular size of the star at the epochs (Sect. |2.2| ). The straight 
lines intersecting at the star's center indicate the biconical 
geometry of the cavities with a projected opening angle Co. 
The cones' axis is tilted by 6 V = 40° with its southern side 
toward the observer, causing the apparently increased cav- 
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Fig. 2. Comparison of the J—H, J—K, and H—K model color images (lower pan els) with o ur latest color images 
of IRC +10 216 (upper panels) obtained in March 2001 with a resolution of 73 mas (Paper III). The observed colors 
were computed for only those pixels where the intensities in H and K are higher than 2 % of the brightest peak, 
whereas the model colors are displayed for the areas brighter that 1 %. The pairs of numbers in the lower corners are 
the maximum and minimum colors. The two dashed circular segments show the radius of the density peak in the main 
dust formation zone of our model at this epoch and the small dashed circles indicate the position and actual size of 
the star. The straight dashed lines intersecting at the star's center outline the conical bipolar cavities of the model. 
Each image shows a 1" x 1" area. North is up and east is to the left. 



ities' opening from uo to uo — 2 arctan (tan (uo/2) cos _1 # v ). 
Two circular segments extended over the angle 7r — uo are 
placed in each image at the density peak position in the 
main dust formation zone of the model. The segments' 
radii increased with time fr om r\ = 28 AU to rs = 33 
AU to rs = 55 AU (Sect. |3.3| ), reflecting the observed 
displacement of component A from the star. 



3.1.1. Motion of the components 

If the bright peaks in our near-IR images of IRC +10 216 
were clumps of gas and dust, their increasing distance from 
the star could be interpreted as real motion of the local 
density enhancements. Our radiative transfer modeling in 
Paper II has shown, however, that the brightest peak A 
seen in all images is not the direct light from the cen- 
tral star. The star is actually located at the position of 
the fainter peak B, whereas component A is the radia- 
tion emitted and scattered in the optically thinner polar 
cavity of the optically thick dust shell. The even fainter 
components C and D were identified with smaller-scale 
deviations of the density distribution of the circumstellar 
environment from axial symmetry. Along the directions to 
these peaks from the central star, circumstellar material 
is less opaque than in other regions of the dense shell. We 



interpreted a difference in the apparent velocities of the 
components (i>a ~ 14 kms -1 vs. vc,r> ~ 5 kins" 1 ) as 
a pure projection effect and deprojected their common 
radial velocity of v r ~ 15 kms -1 for 1995-1997. New 



speckle images presented in [Paper III| revealed a signifi- 
cantly faster motion of component A in 1997-2001, with 
va ~ 21 kms -1 . 

The apparent outward motion can be interpreted ei- 
ther as a real radial expansion of the opaque, dense (spher- 
ical) layer with several optically thinner cavities or as 
a displacement of the dust formation radius due to the 
evaporation of recently formed dust in a hotter environ- 
ment. Increasing temperatures would evaporate grains and 
thus affect the location of the inner dust boundary of 
the envelope, causing its displacement with a velocity 
v T = [<9T(r, t)/dt) [<9T(r, t)/dr}-\ 

Expansion of a dense layer with time-independent pa- 
rameters of dust grains would be inconsistent with the 
star clearly fading since 1997. In fact, any formation and 
expansion of the detached shell would mean that a high 
mass-loss rate has eventually dropped and that the dense 
layer leaves behind lower-density material, as it moves out- 
ward. Since densities and optical depths of such a shell 
would continuously decrease in the process of expansion, 
the latter can only be associated with brightening, not fad- 
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Fig. 3. if -band model intensity profiles compared with 
the intensity cuts (in two orthogonal directions) extracted 
from our images of IRC +10 216 obtained in 1995, 1997, 
and 2001. The unconvolved model intensity cuts through 
the star and the southern cavity show the "true" model 
intensity distributions with a numerical resolution of 5 
mas. 



ing of the star, unless dust grains are still forming in the 
expanding layer, increasing its optical depth. This does 
not seem plausible, because dust formation would quickly 
become inefficient with decreasing gas densities. 

A much better explanation for the decreasing bright- 
ness of the star is to assume a monotonically rising mass- 
loss rate and, hence, higher densities and optical depths 
of the wind in the vicinity of the dust formation radius. 
Note, however, that the increasing M would not form a 
detached dense shell and it is only possible to displace the 
dust boundary if the recently formed dust grains are de- 
stroyed. This would also lead to decreasing optical depths 
of the envelope and brightening of the star, unless dust 
grains continue to form in the wind, making it increas- 
ingly opaque. This scenario is, however, much more likely, 
since the gas density would peak in the dust formation 
zone and it would be increasing with time. 

To distinguish between the real motion of a dense layer 
and the temperature shift of the dust formation zone, we 
need an accurate estimate of the apparent velocities of 
components A, C, and D relative to the star. In fact, one 
can expect an acceleration of the circumstellar material 
by the radiation pressure on dust grains in the neigh- 
borhood of the dust formation zone, not its deceleration. 
If the deprojected radial velocities were not higher than 
the terminal outflow speed Voo « 15 kms -1 observed in 
IRC +10 216, the apparent motion could reflect largely 
an expansion of the detached dense layer in which grains 
are forming. The observed velocities are, however, signifi- 
cantly higher than v^. 

For the assumed distance D = 130 pc, a linear fit to 
the separations of component A ( Paper III ) gives va ~ 18 
kms -1 in the plane of sky, whereas a parabolic fit yields 
va ~ 10 kms -1 in 1996 and suggests 26 kms -1 in 2001 
(va ~ 3 kms -1 yr _1 ). Components C and D moved with 
an average apparent velocity ^c,d ~ 5.5 kms -1 during the 
years of our observations. Interpreting the differences of 
the components' velocitie s in our m odel as pure projection 
effects (cf. Appendix A in Paper II ) and assuming that the 
acceleration of A is real, we can expect for C and D an 
acceleration term vq,b ~ 1 kms -1 yr _1 (^c,d ~ 3 kms -1 
in 1996 and 8 kms -1 in 2001). Within uncer tainties of the 
measured separations, the data presented in Paper III are 
consistent with this picture. For component A, one can 
deproject a radial velocity v r A ~ 19 kms -1 (or v r A ~ 11 
kms -1 in 1996 and 28 kms -1 in 2001, for accelerated 
motion). For components C and D, the velocity is v r c,r> ~ 
17 kms -1 (or v r c,D ~ 9 kms -1 in 1996 and 27 kms -1 in 
2001, if there is acceleration). 

Since the acceleration is most likely real, the expan- 
sion velocity v r ~ 28 kms -1 is now twice higher than the 
terminal wind speed. Taken in context with the increas- 
ing optical depths in the model shell (ryi ~ 30 in 1995, 
rv3 ~ 37 in 1997, and Tv8 ~ 51 in 2001), this strongly sug- 
gests that the observed motions are indeed caused by the 
rapid dust evaporation (vt ~ v r > Voo) due to backwarm- 
ing and higher temperatures in the denser environment 
formed by the increased mass loss (see Sect. pO|). 
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Furthermore, our new calculations demonstrate that it 
may not always be possible to accurately deproject the ob- 
served velocities without a self-consistent radiative trans- 
fer modeling. A reason for this complication is that time- 
dependent optical depth effects arising from the time- 
dependent density distribution (Sect. |3.3| ) play a major 
role in the apparent separations of the observed peaks. 
This is illustrated by the fact that our models, while in 
agreement with the observed distances of component A 
from the star (Figs. |l|, ||), imply significantly higher radial 
velocities of the density peak, v r is ~ 19 kms -1 , v r s8 ~ 27 
kms -1 , and v r i$ ~ 25 kms -1 than simple geometrical 
deprojection would indicate (12 kms -1 , 24 kms -1 , and 
20 kms -1 , respectively). Possible uncertainties in the dis- 
tance to IRC +10 216 cannot decrease the velocities below 

even for the lower limit of 100 pc ( [Becklin et al. 1969| ). 



3.1.2. Near-IR color maps 



Figure Q compares the model and observations in terms 
of the high-resolution J—H, J—K, and H—K color im- 
ages for the most recent epoch of 2001. The speckle color 
images were computed for only those areas of the J, H, 
and K images, where intensities are greater than 2% of 
the peak in the respective bands, whereas the limit for 
model color images was 1%, to display slightly larger ar- 
eas. The obvious difference between the general appear- 
ance of the model and speckle images is mainly due to the 
inhomogeneities of the intensity distributions in the lat- 
ter and a large contribution of components C and D that 
are very prominent at the last epochs, but not included 
in our model. As a consequence, the observed images have 
higher intensity levels over larger areas. One can meaning- 
fully compare only the regions around the bright southern 
cavity at the lates t epochs, in contrast to the H—K image 
of 1997 (P^7T|). 



The color images corroborate the conclusion of 
Paper II that the bright components A, C, and D are 



respectively the cavity and smaller-scale inhomogeneities 
in the dense shell, not dense clumps of dust as it might 
appear from the images alone (Fig. |l|). All the bright com- 
ponents are situated in the blue areas in both the model 
and the observed color images. The cavity A coincides 
with the bluest spot, which is a natural consequence of 
lower optical depths along those directions from the star. 
The bluest, optically thinnest spot is located precisely in- 
side the conical cavity, at the surface of the dense dusty 
shell close to the dust formation radius. Naturally, the star 
is in the red area of the color images, as it was also in the 
previous epochs (Fig. 16 in Paper II). The distribution of 
red and blue color in the observed images would be impos- 
sible to reproduce in any reasonable and accurate model, 
if the star were at the position of the bright southern lobe. 



3.1.3. Shapes of the cavities 

The star is clearly visible in the if -band model images 
at the center of the faint circular halo of dust emission 
from the hottest and densest inner region of the dusty en- 
velope (in 1995 and 1997, Fig. |). The brightest lobe of 
the southern cavity appeares near the edge of the circu- 
lar halo, just at the dust formation zone indicated by two 
dashed circular segments. The intensity level of the halo is 
somewhat higher in the models than in the observed im- 
ages (cf. Fig. ||), causing its prominent appearance. The 
halo is present also in the 1997 K image, where one can 
see the i ncomplete ly circular halo east of the star. As dis- 
cussed in Paper II , the halo is distorted in this image, most 
likely due to density inhomogeneities and patchy circum- 
stellar extinction north and north-west of the star. The 
2001 K image is qualitatively similar to the older image 
obtained in 1997 but strongly distorted by the merging of 
the bright lobe A with component C. Direct stellar light 
has been absorbed in the denser dust shell having much 
larger optical depth (rys ~ 51 vs. Tv3 ~ 36). 

The shape of the cavity has changed between 1995 and 
2001 primarily because of the brightening of component 
C and its merging with A. The latter is somewhat more 
elongated in the model images in the direction orthogonal 
to the symmetry axis (A-B) of the object, which may in- 
dicate that the density distribution or dust properties in 
the cavities in the model differ from the actual parameters 
in IRC +10 216. Despite somewhat different resolutions of 
the images, the size (opening angle) of the outfow cavity 
seems to be larger in 1995 than at later epochs, although 
for 2001 the comparison is complicated by the merger of 
A and C. Our modeling suggests that the cross-section 
of the cavity has indeed decreased with time, when com- 
pared to the cross-section of a cone with fixed opening 
angle. Constrained by the observed if -band speckle im- 
ages, the present models have the opening angle uj of the 
conical cavity decreasing from 46° in 1995 to 36° in 1997 
and to 26° in 2001. We have been unable to find mod- 
els for the first and the last epochs with the same, un- 
changed ou — 36°, that would reproduce intensity ratios of 
the peaks A and B in the observed images. 

This decrease of uj may be interpreted as a kind of "col- 
limation" of the bipolar cavities, although one must bear 
in mind that the "cavities" in reality are not necessarily 
conical, as the model assumes for simplicity. Moreover, the 
details of the intensity distribution in the model images 
are not completely similar to what is shown by observa- 
tions (Fig. ||). The model predicts component A to have 
wider profiles than those of our speckle images, suggesting 
that the density distribution inside the cavities, in contrast 
to the model assumptions, may depend on the polar angle 
or be even more complex. 



3.2. Spectral energy distribution 

Figure compares the model SED with available contin- 



uum fluxes of IRC +10 216 (described in Paper II). Since 
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Fig. 4. Top portions of the model SEDs corresponding 
to the epochs of 1995, 1997, and 2001 com pared to the 
observed fluxes of IRC +10 216 ([Paper Il|) . The beam- 
matched fluxes of the bipolar envelope at all the epochs are 
plotted, as well as the stellar continuum (for 1997). The 
observed broad-band fluxes are shown by different symbols 
for various phases indicated in the legend. If available, er- 
ror bars for the fluxes are drawn, whenever they are larger 
than the symbols. Thick curves (dots) show all available 
spectrophotometry data. The model assumes that we ob- 
serve IRC +10 216 at the viewing angle V = 40° (below 
the midplane). Three arrows at the top abscissa indicate 
the central frequencies of the J, H, and K photometric 
bands. 



the central energy source is a pulsating star, its luminosity, 
effective temperature, and radius are also changing with 
a period of 649 days (Sect. |2.2| ). Our images documented 
the envelope's dynamic evolution over 5.4 years, when the 
material had only time to move by no more than 16 AU 
with the wind velocity of 15 kms -1 . Thus, the observed 
changes have affected only the innermost hot regions of 
the envelope (r ^ 100 AU, T& ^ 300 K) and not influ- 
enced the SED at A ^> 10 /im. For this reason, we display 
in Fig. U only the spectral region 1-50 /im around the 
continuum peak. The entire SED and the prominent dust 
features at 11.3 /im and 27 /im are discussed in detail in 
Paper II . 

The resulting SED displays large phase variations and 
significant non-periodic evolution of the star and the in- 
ner dense envelope of IRC + 10 216. Our model shows that, 
despite the rapid expansion of the dust formation zone, 
the visual optical depth of the dense shell has increased 
from t v1 « 30 to r V3 ~ 37 to r V8 ~ 51 in 1995-2001, 
which implies a significant dust density enhancement in 



the innermost dust formation region (see Sect. 3J3). The 
model predicts that the continuum fluxes at A ^ 5 jam are 
now lower than in 1995 by a factor of ~ 1.5 and that at 
A > 5 /im they are now higher by approximately the same 
factor. The model continuum peak shifted toward longer 
wavelengths from 6.4 jam to 8.6 /im as the increasingly 
dense and opaque dust formation zone expanded radially 
by a factor of 2 (Fig. 0, Sect. The dust emission fea- 
tures at 11.3 /im and 27 /im have now become noticeably 
shallower than several years ago. These model predictions 
can, in principle, be tested by spectrophotometric mea- 
surements. 



3.3. Densities and temperatures 

Figure || displays density and temperature distributions 
in our models of IRC +10 216 for 1995, 1997, and 2001. 
The plots show only the inner regions (r < 10 4 AU) of its 
huge envelope having an outer radius of 6 x 10 5 AU; for 
a discussion of the density and temperat ure structure of 
the entire envelope, we refer to Paper II . The dust tem- 
peratures in Fig. H are plotted for only the smallest dust 
grains of the components considered in our dust model. 
The latter consists of SiC, aggregates composed of sili- 
con carbide and amorphous carbon [SiC,C], and the same 
aggregates covered with t he mantle s of magnesium-iron 
sulfides [SiC,C] + [Mg,Fe]S ( [Paper 13) ). 

The density profile of our previous model was adjusted 
so as to produce the if -band images consistent with the 
observed speckle images of IRC +10 216 (Figs. Q, |). The 
resulting density distribution shown in Fig. || changes with 
time due to the increasing mass-loss rate and forming dust 
in the dense shell expanding away from the central star. 
Between the epochs, the main dust formation zone dis- 
placed from n « 28 AU to r 3 « 33 AU to r 8 « 55 AU. Its 
linear expansion velocity in our model, v r ss ~ 27 kms -1 , 
is almost twice larger than the general outflow speed of 
the envelope (Sect. |3.1.l| ). 

The total density distributions shown in Fig. |H| assume 
that the dust-to-gas mass ratios p&j p of different dust 
components in IRC +10 216 do not change in time and 
that they are equal to the values adopted in our p revi- 
ous time-independent model (see Table 2 in Paper II). An 
inspection of Fig. || shows that the densities in the dust 
formation zones of different dust materials became higher 
with time, implying that the mass-loss rate M increased, 
too. For an illustration, we plotted p oc r -2 density pro- 



files through the dust and total density peaks of 1995. 
These lines would trace the displacement of the peaks with 
time, if M, u, and pd/ p remained constant. Contrasting 
with the expected density distribution of such an outflow, 
the model densities rise significantly above the pure ex- 
pansion line. This can be interpreted in terms of the in- 
creasing mass-loss rate or of the increasing degree of dust 
condensation or both. 

In the standard picture of mass-losing AGB stars, the 
radiation pressure acceleration of dust and gas takes place 
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Fig. 5. Model densities (left panel) and dust temperatures (right panel) in the inner envelope of IRC +10 216. Dust 
densities and temperatures are displayed for the smallest grains only (a = 0.01 pm) of all dust components. Both total 
(gas+dust) density and dust density distributions in the midplane and in bipolar cavities along the symmetry axis are 
shown here. The total density profiles were plotted under the assumption that pd/ P of different dust components do 
not change with time. Three vertical dotted lines in the middle of the plot mark the amorphous carbon condensation 
zones in 1995, 1997, and 2001. Jumps in densities and temperatures in the area are due to sequential formation of 
new dust components (Paper II). Dash-double-dotted lines show a p cx r~ 2 density profile that would form in the case 



of constant M and v. The additional labeling above the lower abscissa shows expansion times assuming a constant 
outflow velocity v = 15 kins" 1 . The two vertical lines drawn on the left side of the diagrams (labeled with an asterisk) 
indicate stellar radii at the luminosity phases <p = 0.0 and <p = 0.5. 



in the dust formation zone. In reality, different dust com- 
ponents form within a considerable range of radial dis- 
tances from the star (cf. |Paper I]| ), although the general 
picture of the acceleration of the newly-formed dust grains 
due to the radiation pressure remains valid. The velocity 
v r 38 « 27 km s _1 of the wind material close to the dust for- 
mation zone, which is higher than the generally accepted 
value of Voq ^15 kms -1 , does not fit readily into this con- 
text, as it would imply significant deceleration of the wind 
material. Although such a deceleration may well be possi- 
ble in a compl ex, time-dependent outflow hydrodynamics, 
m Sect. |3.1.1| we have interpreted this fast expansion as 
not the real outflow motion, but rather a displacement of 
the dust formation radius due to the backwarming effect 
in the optically thicker environment. 

In fact, the radial dust temperature distributions of 
our model (Fig. ||) demonstrate a very significant heating 
of the grains within the dense shell. The model density 
profile just outside the carbon nucleation zone steepens 
with time from pdi cx r -3,8 to pd3 cx r -4,5 to pas cx r -8 
between the three epochs, causing higher values and much 
steeper distributions of dust temperatures in that region. 
When they jump over the condensation temperature of the 
dust material, evaporation of the dust grains is a natural 



and inevitable consequence. Dust sublimation can quickly 
shift the dust formation zone outward, thus bringing the 
increased densities in balance with the condensation tem- 
peratures of carbon dust, Tq ~ 10 3 K, however, at a larger 
radius. Such displacement of the dust density peak due to 
the self-regulating processes of dust formation, heating, 
and destruction can easily be mistakenly interpreted as a 
radial expansion of the wind material. 

3.4. Mass-loss rate and dust formation 

The time-dependent mass- loss rate M of IRC +10 216 can 
readily be derived from the density distribution of our 
models, assuming that the outflow velocity v = 15 kms -1 
is constant across the envelope and that pd/ p in the en- 
velope does not change with time (Fig. |6|). Note that the 
radiative transfer models like this predict only dust den- 
sities and convert them into gas densities by adopting a 
specific value (usually arbitrary) of the dust-to-gas mass 
ratio. However, pd/ p is not a free parameter in our model 
of IRC +10 216 (see [Paper lip . 

If we assume that pd/ p remained constant during the 
envelope's expansion, then the rising densities displayed 
in Fig. |5| would imply monotonically increasing mass-loss 
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rates illustrated in Fig. g. In this case, our model predicts 
that the mass-loss rate has increased by a factor of 7 dur- 
ing the 5.4 years of our imaging, from Mi ~ 8 x 10 -5 
M yr- 1 in 1995 to M 3 « 1.4 x 1(T 4 M yr" 1 in 1997 to 
M 8 ~ 5.8 x 10" 4 Mqyt- 1 in 2001. Such changes would 
amount to a rate of ~ 10 -4 Mq yr -2 . Assuming that the 
latter is constant and that it takes 20 years for the ma- 
terial leaving the star to reach the zone of carbon dust 
formation (at ~ 55 AU, Fig. ||), we would conclude 
that the mass-loss rate from the stellar surface may have 
been recently as high as 2.6 x 10 -3 Mq yr -1 . 

If we assume instead that M remained constant during 
the envelope's expansion, then the total density profiles 
must differ from those plotted in Fig. |[ The total den- 
sity peak of 1995 should evolve then along the p oc r -2 
track. This would imply that pd/ P has increased by the 
same factor of 7, from pdi/p ~ 9-4 x 10 -4 in 1995 to 
Pds/p ~ 1-65 x 10~ 3 in 1997 to p d8 /p « 6.6 x 10~ 3 in 
2001, i.e. that the new dust formation continues as the 
dense gas flows outward. This seems to be natural, since 
the dust forma tion process takes considerable time (see, 
e.g., Paper III ). In this case, the mass- loss rate would re- 
main at a lower value of 8 x 10 -5 Mq yr -1 but pd/ p would 
be rather high at present, approaching 0.7%. 

It seems very likely, however, that in IRC +10 216 both 
M and pd/ P are increasing with time. If this indeed is the 
case, their values can be not as high as in the above two 
scenarios. For example, if their contributions were equal, 
then we would derive ~ 3.3 x 10 -4 MQyr -1 and 
Pds/p * 0.004. 

4. Conclusions 

We presented results of two-dimensional radiative transfer 
modeling of IRC +10 216 at selected moments of its recent 
evolution, which correspond to three epochs of our series 
of high-resolution spec kle imag e s of the o bject recorded 
over the last 5.4 years ( Paper I ; Paper III ). The imaging 
revealed dynamic evolution of the subarcsecond dusty en- 
vironment of IRC +10216 on a very short time scale, of 
the order of one year. The goal of the present study was to 
t est predi ctions of our recent time-independent modeling 
( Paper II ) which derived the structure and physical prop- 
erties of IRC +10 216 at a single epoch (January 1997). 
Now that the entire sequence of 8 high-resolution near- 
IR images has documented complex changes in the inner 
bipolar dusty shell of the carbon star, we took these tem- 
poral constraints into account to see, whether the static 
model is consistent with the fast evolution seen in our im- 
ages. The new modeling allowed us to make a quantitative 
physical interpretation of what goes on in IRC +10 216. 

Our previous static model of IRC +10 216 ( Paper II ) is 
consistent with the new constraints, provided that small 
modifications are made to its density profile in the dust 
formation zone and to the opening angle of the bipolar 
cavity. The new modeling has shown that the cavity has 
been shrinking from uj\ « 36° to us ~ 26° and the den- 
sity distribution across it has been changing during the 
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Fig. 6. Mass- loss history of the inner envelope of IRC 
+10 216 for October 1995, January 1997, and March 2001, 
derived from the total density distributions displayed in 
Fig. ||. The total mass-loss rates in the midplane and along 
the model symmetry axis (in polar cavities) are plotted, 
assuming that the outflow velocity v = 15 kms -1 is con- 
stant throughout the envelope and that pd/ p of all dust 
components does not change with time. The dash-double- 
dotted line corresponds to a constant M spherical outflow 
resulting from the p oc r~ 2 density profile (Fig. ||). 



5.4 years of our imaging. We are witnessing a dynamic 
episode of rising mass loss from the central star, from 
M « 10 -5 Mq yr -1 to the present value of M « 3 x 10 -4 
Mq yr -1 , which probably started ~ 50 years ago. If the 
current rate of the increase of M is constant, the mass- 
loss rate from the stellar surface may be now as high as 
2.6 x 10~ 3 M^yr- 1 . 

A compact dense shell with bipolar cavities has formed 
around the star as a result of the rapid increase of the mass 
loss by IRC +10 216, causing the observed rapid evolution 
documented in our images. The higher mass loss produces 
favorable conditions for dust formation in the increasingly 
dense inner envelope expanding outward with the outflow 
velocity v ~ 15 kms -1 . Larger amounts of dust increased 
the optical depths, obscuring the central star, whereas 
the optically thinner cavity remained relatively unaffected. 
Due to backwarming, temperatures in the dust formation 
zone became higher, shifting the latter to larger radii at 
the velocity vt ~ 27 kms -1 and mimicking the real out- 
flow motion of the circumstellar shell material with that 
speed. One can predict that the star will remain obscured 
until M starts to drop back to lower values in the dust 
formation zone. Within a few years from that moment, 
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we could be witnessing the star reappearing, whereas the 
cavities becoming relatively fainter. 
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